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ABSTRACT 

MA% is a new tool to estimate parameters from stellar spectra. It is based on the maximum likelihood method, with the likelihood 
compressed in a way that the information stored in the spectral fluxes is conserved. The compressed data are given by the size of the 
number of parameters, rather than by the number of flux points. The optimum speed-up reached by the compression is the ratio of the 
data set to the number of parameters. The method has been tested on a sample of low-resolution spectra from the Sloan Extension for 
Galactic Understanding and Exploration (SEGUE) survey for the estimate of metallicity, effective temperature and surface gravity, 
with accuracies of 0.24 dex, 130K and 0.5 dex, respectively. Our stellar parameters and those recovered by the SEGUE Stellar 
Parameter Pipeline agree reasonably well. A small sample of high-resolution VLT-UVES spectra were also used to test the method 
and the results have been compared to a more classical approach. The speed and multi-resolution capability of MAx combined with 
its performance compared with other methods indicates that it will be a useful tool for the analysis of upcoming spectral surveys. 
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1. Introduction 

The astronomical community has conducted many massive 
surveys of the Universe, and many more are either ongoing or 
planne d. Recently comp leted is the Sloan Digital Sky Survey 
(SDSS. lYorket.al.l2 0001. notable for assembling in a consistent 
manner the positions, photometry and spectra of millions of 
galaxies, quasars and stars, and covering a significant fraction 
of the sky. The size of such a survey allows us to answer many 
questions about the structure and evolution of our universe. 
More locally, massive surveys of stars have been undertaken 
to reveal their properties, for example t he Geneva-Copenhagen 
Survey for the solar neighbor hood ([Nordstrom et al.l 2004) 
and the ELODIE library dPrugniel et alj|2007l) . These surveys 
provide a copious amount of stellar photometry and spectra 
from the solar neighborhood, broadening the range of stellar 
types studied. The SEGUE catalog - Sloan Extension for 
Galactic Understanding and Exploration ( Yanny et al. 2003 - a 
component of SDSS, contains additional imaging data at lower 
galactic latitudes, to better explore the interface between the 
disk and halo population. 

By statistically analyzing the properties of these stars such 
as chemical abundances, velocities, distances, etc, it has been 
possible to match the structure and evolution of the Milky Way 
to the current generation of galaxy formation models ( Juric et all 
l20Q8t llvezic et al.ll2008t iBond et al.ll2009l IWvsell2006l) . " 

Studies of large samples of stars of our galaxy are crucial 
tests for the theory of the structure and formation of spiral galax- 
ies. Even though the statistics given by SEGUE or the Geneva- 
Copenhagen survey agree with the models and simulations of 
galaxy formation, there is some contention over whether the ac- 
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curacies of measurements ade quately support the co nclusions. 
New surveys such as RAVE dSteinmetz et al.l [2006) and Gaia 
dPerrvman et alj I2001I) will give ex tremely high accu racies in 
velocities and positions, LAMOST (IZhao et alJl2006|) and part 
of the SDSS-III project, APOGEE and SEGUE-2 iRockosi et al.l 
2009), will give high-quality spectra and therefore stellar param- 
eters and chemical abundances of millions of stars over all the 
sky. These more accurate properties and larger samples will in- 
crease our knowledge of the Milky Way and answer wider ques- 
tions about the formation and evolution of spiral galaxies. 
As data sets grow it becomes of prime importance to create effi- 
cient and automatic tools capable of producing robust results in 
a timely manner. A standard technique for estimating parameters 
from data is the maximum-likelihood method. In a survey such 
as SEGUE, which contains more than 240,000 stellar spectra, 
each with more than 3000 flux measurements, it becomes ex- 
tremely time-consuming to do spectral fitting with a brute-force 
search on the multi-dimensional likelihood surface, even more 
so if one wants to explore the errors on the recovered parame- 
ters. 

Efforts to reduce the computational burden of characterizing the 
likelihood surface include Markov-Chain Monte-Carlo methods 
(Ball ot etalJl2 006). where a chain is allowed to explore the like- 
lihood surface to determine the global solution. An alternative 
method is to start at an initial estimate and hope that the likeli- 
hood surface is smooth enough that a gradient search will con- 
ve rge on the solution. An exam ple of t his approach can b e found 
in lAllende Prieto et. all (I2006I) and iGrav etal l (l200lh. based 
on th e Nelder-Mead downhill simplex method dNelder & Mead 
1965), where the derivatives of the likelihood function give the 
direction toward the maximum. In both cases, the likelihood 
is evaluated only partially, leaving large amounts of parameter 
space untouched. Moreover, the time needed to find the maxi- 
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mum depends on the starting point and the steps used to evaluate 
the next point. 

Another parame ter estimation meth od u tilizes neural net- 
works , for example ISnider et al.l (12001 1) and iRe Fiorentin et al.1 
(l2007h . While these non-linear regression models can obtain 
quick accurate results, they are entirely dependent on the quality 
of the training data, in which many parameters must be previ- 
ously known. Careful attention must also be paid to the sampling 
of the model grid which is used to generate the neural network. 

An alternative approach to spectral analysis is the Massively 
Optimized Parameter Estimation and Data compression method 
(MOPED Q iHeavens etail 12000b . This novel approach to the 
maximum likelihood problem involves compressing both data 
and models to allow very rapid evaluation of a set of parame- 
ters. The evaluation is fast enough to do a complete search of the 
parameter space on a finely resolved grid of parameters. Using 
carefully constructed linear combinations, the data are weighted 
and the size is reduced from a full spectrum to only one num- 
ber per parameter. This number, with certain caveats (discussed 
later), contains all the information about the parameter con- 
tained in the full data. The meth od has been succe ssfully applied 
in the fields of CMB analysis dBond et al.l 1998), medi cal im - 
age registration and galaxy shape fitting dToieiro et al.l 12007). 
A complete background of the development of the MOPED al- 
gorithm can be fou nd in Tegmark et all (Il997|. hereafter T 97), 
IHeavens et al] d2000l hereafter H00) and lPanter etall d2003l) . 

We present a new derivative of MOPED, MA^ (MAssive 
compression of x 2 ), to analyze stellar spectra. To estimate the 
metallicity history of a galaxy for instance, MOPED needs mod- 
els where the spectra are the sum of single stellar populations. 
Here we study one population, meaning the metallicity has to 
be estimated from the spectrum of each single star. It is there- 
fore necessary to develop a specific tool for this task: MA^. In 
Sect. |2]we describe the method, giving a summary of the deriva- 
tion of the algorithm. We then test the method in Sect. [3] on 
the basic stellar atmosphere parameter estimation of a sample 
of low-resolution stars of SEGUE, where we compare our re- 
sults with the SEGUE Stellar Parameter Pipeline in Sect. [4] and 
finally check the method for a small sample of high-resolution 
spectra in Sect. [5] where we compare our results with a more 
classical approach for a small sample of VLT-UVES spectra. A 
summary of our conclusions is given in Sect. [6] 

2. Method 

In this section we describe our method to treat the likelihood sur- 
face. We first review the classical maximum-likelihood method, 
where a parametric model is used to describe a set of data.Then 
we present the algorithm that compresses the data and hence 
speeds the likelihood estimation. Finally we show the proof for 
the lossless nature of the compression procedure. 



2. 1 . Maximum-likelihood description 

Suppose the data (e.g. the flux of a spectrum) are represented 
by N real numbers x\,x>i, ...,Xn, which are arranged in an 
/V-dimensional vector X. They are the flux measurements at N 
wavelength points. Assume that each data point x, has a signal 
part yu,- and a noise contribution <x, 

X i =/l i + (T i . (1) 



1 MOPED is protected by US Patent 6,433,710, owned by The 
University Court of the University of Edinburgh (GB) 

2 http://www.blackfordanalysis.com 



If the noise has zero mean, (X) = (//), we can think of X as 
a random variable of a probability distribution L(X, 0), which 
depends in a known way on a vector of m model parameters. 
In the present case m = 3 with 



= (8 l ,0 2 ,e 3 ) = ([Fe/H], T eff , log g), 



(2) 



with the iron abundance [Fe/H], effective temperature T e ff and 
surface gravity logg. A certain combination of the single param- 
eters 6j (j = 1,2,3) also produces a theoretical model F = F(0), 
and we can build an ra-dimensional grid of theoretical models 
by varying the 8j parameters. If the noise is Gaussian and the 
parameters have uniform priors, then the likelihood 
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gives the probability for the parameters, where cr 2 ^ is the aver- 
aged square of the noise of each data point and 
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The position of its maximum estimates the set of parameters 0o 
that best describe the data X. The fit between F(0o) and X is 
good if x 1 < DOF, where DOF = N — m are the degrees of free- 
dom. In the most basic form one finds the maximum in the likeli- 
hood surface by exploring all points in the parameter space, with 
each likelihood estimation calculated using all N data points. 
This procedure is of course very time-consuming if N and m are 
large. 

2.2. Data compression 

In practice not all data points give information about the param- 
eters, because either they are noisy or not sensitive to the pa- 
rameter under study. The MOPEDalgorithm uses this knowledge 
to construct weighting vectors which neglect some data without 
losing information. A way to do this is by forming linear combi- 
nations of the data. 

Let us compress the information of a given parameter 9j from 
our spectrum. The idea is to capture as much information as pos- 
sible about this particular parameter. We define our weighting 
vector bj = (fe/,i,fc/,2, bjjf) as 



1 1 <9F 



(5) 



where |bj| = Jbj ■ bj. This definition assures that each data 
point of the spectrum is less weighted if <x is large and more 
weighted if the sensitivity in the flux - meaning the derivative 
of F respect to the parameter Oj - is high. To create this vector 
we need, on one hand, the information about the behavior of the 
parameters from a theoretical fiducial model, and on the other 
hand, the error measurement from the data. 

Figure Q] shows the example of the weighting vectors for 
the parameters of the set (f2]) using a fiducial model with 
[Fe/H] = -2.0, T eff = 6000 K and logg = 4.0. The upper panel 
corresponds to the synthetic spectrum at the SDSS resolving 
power (R = 2000) and the panels B, C and D to the weighting 
vector for metallicity, effective temperature and surface gravity, 
respectively. This figure is a graphical representation of Eq. (0, 
and helps to better visualize the weight in relevant regions of the 
spectrum. 
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Fig. 1. (A) Fiducial model with parameters [Fe/H] = -2.0, 
T e ff = 6000 K and \ogg = 4.0. Other panels indicate the weight- 
ing vectors according to Eq. (|5]l for metallicity (B), temperature 
(C) and surface gravity (D). 



For metallicity (panel B), the major weight is concentrated 
in the two peaks at 3933 and 3962 A, corresponding to the Call 
K and H lines seen in the synthetic spectrum in panel A. The 
second important region with weights is close to 5180 A, corre- 
sponding to the Mglb triplet. A minor peak is seen at 4300 A, 
which is the G-band of the CH molecule. These three features 
are metallic lines, therefore show a larger dependency on metal- 
licity than the rest of the vector, which is dominated by a zero 
main weight. 

For the temperature (panel C) the greatest weight is focused 
on the hydrogen Balmer lines at 4101, 4340 and 4861 A. Peaks 
at the previous metallic lines are also seen, but with a minor 
amplitude. Finally, the surface gravity (panel D) presents the 
strongest dependence on the Mgl triplet. The wings of the Call 
and Balmer lines also present small dependence. As in the case 
for metallicity, the rest of the continuum is weighted by a mean 
close to zero. 

The weighting is done by multiplying these vectors with the 
spectrum, where the information kept from the spectrum is given 
by the peaks of the weighting vector. We define this procedure 
as "compression", where the information about the parameter 0j 
is expressed as 
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Fig. 2. Correlations between the Fisher matrix values obtained 
from the full and compressed likelihoods for metallicity (upper 
panel) and effective temperature (lower panel) for 75 randomly 
selected stars of SEGUE. The line corresponds to the one-to-one 
relation. 



with X the spectrum. With the weighting vector defined in 
Eq. ([5]) we can compress data (X) and model (F) by using Eq. ©. 
Because bj and X (or F) are A^-dimensional, the product yj in 
(0 is a number, which stores the information about the parame- 
ter 

In order to perform a compression for all the parameters si- 
multaneously, we require that y* is uncorrected with yj, with 
j + k. This means that the b-vectors must be orthogonal, i.e. 
bt ■ bj = 0. Following the procedure of H00 we find the other 
yk — bk ■ X by the Gram-Schmidt ortogonalization 



bk = in 

|b*| 



9=1 



q] b q 



(7) 



For m — 3 we have the numbers yi_x, yi.x, y3,x, f° r the data and 
yi,F>y2,F,y3,F for the models, corresponding to the parameters of 
given in Eq. (f2]l 

Goodness of fit 

The definition of y 1 , which is the sum of the differences of the 
data and the model, motivates us to define our goodness of fit. 
The "compressed" xl is me sum of the differences of the com- 
pressed data and model 



yj* = ^ • X, 



j - 1,2,3 



(6) 
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-(yjx-yjf) 
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which gives the probability for the parameter Gj. 
1 



Lcj(x,ej) 



(2n) N ' 2 Vo-; 



exp ■ 



(9) 



Because the yj numbers are by construction uncorrected, the 
"compressed" likelihood of the parameters is obtained by mul- 
tiplication of the likelihood of each single parameter L c j(X, Gj) 



L c (X,&) = Y\L cJ (K,ej) 



1 



((2n) N / 2 



where 



1 m 

xl = — Yj(y.j,x-yj,F) 2 



exp^, (10) 



(11) 



is the compressed^ 2 . 

As in the classical approach, the peak of the compressed 
likelihood estimates the parameters 0q that generate the model 
F(0o), which reproduces best the data X. The advantage of using 
the compressed likelihood is that it is fast, because the calcula- 
tion of xl needs only m iterations and not N as in the usual y 2 
computation. The search of the maximum point in the likelihood 
is therefore done in an m-dimensional grid of y-numbers and 
not of iV-length fluxes. 

2.3. No loss of information by the compression 

The Fisher Information Matrix describes the behavior of the like- 
lihood close to the maximum. Here we use it only to show the 
lossless compression offered by the MOPED method, but a more 
extensive study is given in T97 and H00. 

To understand how well the compressed data constrain a solution 
we consider the behavior of the logarithm of the likelihood L 
near the peak. Bellow we denote a generic L that can correspond 
to the compress likelihood L c or the full one Lf, because the 
procedure is the same. In the Taylor expansion the first deriva- 
tives d(ln L)/d9j vanish at the peak and the behavior around it is 
dominated by the quadratic terms 



InL = InL, 



1 d 2 lnLi 



peak 
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The likelihood function is approximately a Gaussian near the 
peak and the second derivatives, which are the components of 
the inverse of the Fisher Matrix Tjk, measure the curvature of 
the peak. Because the dependence of the parameters is not cor- 
related, the matrix is diagonal with values 



5 2 lnL 
86) ■ 



(13) 



In Fig.|2]we plotted the values !F[Fe/H] and Tj CIT , correspond- 
ing to the parameters of metallicity and effective temperature, 
respectively, for likelihoods of 75 randomly selected stars of 
SEGUE, using the classical (full) and the compressed (comp) 
dataset. The line corresponds to a one-to-one relation. The val- 
ues correlate well, as expected. 

Figure[3]shows full and compressed likelihood surfaces of metal- 
licity and temperature for four randomly selected stars, where 
eight equally spaced contour levels have been plotted in each 
case. This is another way to visualize the statement that the 
Fisher Matrices correlate. The curvature of both likelihoods is 
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Fig. 3. Likelihoods of the compress data set (left) and full data 
set (right) in the parameter space of effective temperature and 
metallicity for four randomly selected SEGUE stars. In each 
panel eight equally spaced contour levels are plotted. The trian- 
gle and diamond correspond to the maximum point of the com- 
pressed and the full data set, respectively. Both maxima lie in the 
first confidence contour level, meaning that both data sets reach 
the same solution for these two parameters. 

the same close to the peak, as predicted. In Fig. [3] we also over- 
plotted the maximum point of the compressed likelihood as a 
triangle and the maximum point of the full one as a diamond. It 
can be seen that the maximum points of the compressed and full 
data set lie inside the first contour level of the likelihoods. This 
shows that the compressed data set gives the same solution as 
the full one - even after the dramatic compression. 

In the bellow parts the numerical values of x 1 correspond to 
the classical reduced^ 2 . 



3. Implementation on low-resolution spectra 

The sample of stars used to test the method are F-G dwarf stars. 
Based on the metallicity given by the SEGUE Stellar Parameter 
Pipeline, we chose metal-poor stars as will be explained below. 
This choice allows us to avoid the problem of the saturation of 
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metallic lines such as Call K (Beers et al . 1999), which is a very 
strong spectral feature that serves as a metallicity indicator in the 
low-resolution spectra from SEGUE. These stars fall in the tem- 
perature range where B aimer lines are sensitive to temperatur e 
and the spectral lines are not affected by molecules dGravll 19921) . 
With these considerations it is correct to assume that the spec- 
tra will behave similar under changes of metallicity, temperature 
and gravity. This allows the choice of a random fiducial model 
from our grid of models for the creation of the weighting vec- 
tors, which will represent well the dependence of the parameters 
in all the stars. 

3.1. Data 

We used a sample of SEGUE stellar sp ectra (lYannv et al.ll2 009). 
part of the Seventh Data Release (DR7, Abazaiian et alJ2009l) of 
SDSS. The survey was performed with the 2.5-meter telescope 
at the Apache Point Observatory in southern New Mexico and 
contains spe ctra of ~ 240, 000 sta rs. All the spectra have ugriz 
photometry dFukugita et al.lll9 96). From a color-color diagram 
we iden tified the dwarf sta rs and selected them using the sql- 
Loader (Thakar et al. 2004) under the constraints that the object 
must be a star and have colors of 0.65 < u — g < 1.15 and 
0.05 < g - r < 0.55. We verified the spectral type by also down- 
loading the values 'seguetargetclass' and 'hammerstype', which 
classify our stars mainly as F-G dwarfs. We also constrained the 
metallicity to be in the range of -999 < [Fe/H] < -0.5 . The 
latter valu es were taken fro m the new SEGUE Stellar Par ameter 
Pipeline dLee et al.1 12008allb1; lAllende Prieto et al.1 120081) . were 
the -999 indicates that the pipeline has not estimated the metal- 
licity of this particular star. 

Our final sample contains spectra of 17,274 stars with a resolv- 
ing power of R 2000 and a signal-to-noise up to 10. The wave- 
length range is 3800 - 9200A and the spectra are with absolute 
flux. 

3.2. Grid of models 

The synthetic spectra were created with the synthesis code 
SPECTRUM dGrav & Corballvll 19941) . which uses the new fully 
blanketed stellar atmosphere models of iKuruczl (fl 992) and com- 
putes the emergent stellar spectrum under the assumption of lo- 
cal thermal equilibrium (LTE). The stellar atmosphere m odels 
assume the solar abundances of lGrevesse & Sauvall ([T998) and a 
plane parallel line-blanketed model structure in one dimension. 
For the creation of the synthetic spectra, we set a microturbu- 
lence velocity of 2 km/s, based on the atmosphere model value. 
The line-list file and atomic data were taken directly from the 
SPECTRUM webpageQ In these files, the lines were taken from 
the NIST Atomic Spectra Database Q and the Kurucz webside0. 
No molecular opacity was considered in the model generation. 

We created an initial three dimensional grid of synthetic 
spectra st arting from the ATLAS9 Grid of stellar atmosphere 
models of ICasteili & Kuruczl (l2003h by varying the parameters 
of Eq. They cover a wavelength range from 3800 to 
7000 A in steps of AA — I A, based on the wavelength range of 
SDSS spectra together with our line-list. This wavelength range 
is broad enough for a proper continuum subtraction, as described 
in Sect. 13.31 The spectra have an absolute flux and were finally 



3 http://wwwl.appstate.edu/dept/physics/spectrum/spectrum.html 

4 http://physics.nist.gov/PhysRefData/ASD/index.html 

5 http://kurucz.harvard.edu/linelists.html 
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Fig. 4. Radial velocities found using minimum fluxes of strong 
lines given by Table Q] (MAX) compared with those of the 
SEGUE database. The difference of the radial velocities given 
by the SEGUE database from those obtained by our method is 
indicated in the legend as offset, with its standard deviation as 
scatter. 



smoothed to a resolving power of R — 2000, according to SDSS 
resolution. 

In order to have a finer grid of models, we linearly in- 
terpolated the fluxes created for the initial grid. It has mod- 
els with -2.5 < [Fe/H] < -0.5 with A[Fe/H] = 0.1 dex, 
5000 < T eff < 8000 with AT eff = 50 K and 3.5 < logg < 5 
with Alogg = 0.1 dex. The scaling factor for the normalization 
varies linearly from 0.85 to 1.15 in steps of 0.01 (see below). 
Our final grid has 21 x 61 x 16 x 30 = 614,880 models. We 
set a linearly varying [a/Fe] of [a/Fe] = [0.2,0.4] for stars in 
the metallicity range of -0.5 > [Fe/H] > -1.0 and [ar/Fe] = 0.4 
for stars with [Fe/H] < -1.0, to reproduce the a bundance of a 
elements in the Milky Way, as in lLee et aL (2008a). These vary- 
ing [a/Fe] abundances were also calculated with interpolation of 
fluxes created from solar and a- enhanced ([a/Fe] = 0.4) stellar 
atmosphere models. We did not include the a abundances as a 
free parameter, because we aim to estimate the metallicity from 
lines of a elements (Ca, Mg), meaning we would obtain a de- 
generacy in both parameters. We prefered to create a grid where 
the a abundances were already 'known'. 

The grid steps between the parameters are smaller than the 
accuracies expected from the low resolution SEGUE data. The 
extra time required to calculate the compressed likelihood in this 
finer grid is not significant, and retaining the larger size allows 
us to demonstrate the suitability of the method for future more 
accurate data capable of using the full grid. 

3.3. Matching models to data 

Data and model needed to be prepared for the analysis. First of 
all, the data needed to be corrected from the vacuum wavelength 
frame of the observations to the air wavelength frame of the lab- 
oratory. Secondly, we needed to correct for the Doppler effect zc. 
This was done by using the flux minimum of the lines indicated 
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Table 1. Strongest lines in F-G dwarf stars. The first column 
indicates the name of the line and the second its wavelength po- 
sition in A. 



line 


A A 


Ca II K 


3933 


C&IIH 


3962 


m 


4101 


Hy 


4340 


W 


4861 


Mg b triplet 


5183 



in Table Q] except for the Mglb triplet, because this feature is 
not seen clearly enough in every spectrum, driving in our auto- 
matic zc calculation to unrealistic values. A comparison of the 
radial velocity found with this method with the value given by 
the SEGUE database f is shown in Fig. [4] The difference of the 
radial velocities (SEGUE - MAX) has a mean (offset) of -20.15 
km/s with standard deviation (scatter) of 7.85 km/s. The effect of 
this difference on the final parameter estimation is discussed in 
Sect|4] Once our data were corrected for the Doppler effect, they 
had to be interpolated to the wavelength points of the models to 
get the same data points. 

For automatic fitting of an extensive sample of stars with a large 
grid of models, we decided to normalize the flux to a fitted con- 
tinuum. In this way the dependence of the parameters is this is 
not past, it is a fact that happens now too) concentrated mainly 
on the line profiles. The choice of the normalization method is 
a difficult task, because none of them is perfect. It becomes es- 
pecially difficult in regions where the spectrum shows too many 
strong absorption features, as the Call lines of TableQ] Using the 
same method for models and data, difficult regions behave sim- 
ilar in both cases, which allowed us a final fitting in these com- 
pl icated regions as well. We a dopted the normalization method 
of lAllende Prieto et. al.l d2006l) . because it works well for the ex- 
tended spectral range of SDSS spectra. It is based on an itera- 
tively polynomial fitting of the pseudo-continuum, considering 
only the points that lie inside the range of 4cr above and 0.5cr 
below. Then we divided the absolute flux by the final pseudo- 
continumm. Due to noise, the continuum of the final flux may 
not necessary be at unity. We included a new degree of freedom 
in the analysis that scaled this subtracted continuum. This scal- 
ing factor for the normalization is then another free parameter, 
which means that we have finally four parameters to estimate - 
three stellar atmosphere parameters and the scaling factor. The 
final step is to choose the spectral range for the analysis. 

3.4. Compression procedure 

We chose the fiducial model with parameters [Fe/H] = 
-2.0,T eff = 6000 K, logg = 4.0 and a s = 1, to calculate 
the weighting vectors, b[F e /H], bx^, bi og? and b aj , corresponding 
to metallicity, effective temperature, gravity and scaling factor 
for normalization, respectively. Then, we calculated the set of 
j-numbers using Eq. © for each point in the grid of synthetic 
spectra by projecting the b- vectors onto the spectrum, resulting 
in a four dimensional y-grid, with every point a single number. 
With the respective b-vectors we calculated the y-numbers for 
the observed spectra, the expression (TTTT) for every point in the 
grid and finally we found the minimum value in the grid which 
corresponds to the maximum point of the compressed likelihood. 

To refine our solution we found the "real" minimum^ 2 with 
a quadratic interpolation. For the errors, we looked at the models 
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Fig. 5. Example of the fit between a randomly selected SEGUE 
star (crosses) and a synthetic spectrum (red dashed line). The 
legend indicates the stellar atmosphere parameters of the model 
and the value of the reduced^ 2 of the fit. The upper panel is the 
fit using all the points of the spectral region [3850,5200] A. The 
lower panel is the fit using only the data points where the lines 
of Table[TJare located. 



within the confidence contour of x = X m i„ + A* w ith At = 
4.72 corresponding to the l<x error in a likelihoo d with four free 
parameters (for further details see chapter 14 of Press 1993b. 



4. Application to SEGUE spectra 

We chose to use the wavelength range of [3850,5200] Aas in- 
formation about metallicity, temperature and gravity, which is 
available in the lines listed in Table [TJ These spectral lines are 
strong in F-G stars, meaning they can be identified at low reso- 
lution without difficulty. The wings of Balmer lines are sensitive 
to tempera ture and the win gs of strong Mg lines are sensitive 
to gravity dFuhrmannlll998h . Beacuse iron lines are not strong 
enough to be distinguished from the noise of our spectr a, the Call 
K and Mglb fines are our indicators of metallicity dBeers et al] 
120001: lAllende Prieto et. afteOOdlLee et alJl2008ah . 

It is important to discuss the carbon feature known as the G- 
band at 4304 A in our spectral window. The fraction of carbon- 
enhanced metal-poor ( CEMP) stars is exp ected to increase with 
decreasing metallicity (Beer s et ail 1 199 2) to about 20% at be- 
low [Fe/H] = -2.0 dLucatello et al.l l2006).possiblv implying a 
stron g carbon feature in the observed spectrum. As discussed 
by M arsteller et alj (2009), a strong G-band may also affect the 
measurement of the continuum at the Call lines, which could 
result in an underestimation of the stellar metallicity. We com- 
mented in Sect. [2] on the influence of the lines in the weighting 
vectors used for the compression and we saw that the G-band 
also displayes minor peaks (see Fig. [TJ). The role of this minor 
dependence compared with those from the Call , Balmer and 
Mgl lines can be studied by comparing the final parameter es- 
timation when using the whole spectral range or only those re- 
gions with the lines of Table QJ where the G-band is not con- 
sidered. A further motivation for performing an analysis of an 
entire range against spectral windows is also explained below. 
The implications of this tests in terms of fraction of CEMP stars 
is discussed below. 
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Fig. 6. Upper panel: Metallicity (left), effective temperature (middle) and surface gravity (right) obtained using MA% for the entire 
spectral range of [3900, 5200] A (whole) compared to that with selected spectral windows (win) for a sample of 150 randomly 
selected stars. The offset (mean difference of "win - whole") of the results and its scatter (standard deviation) is indicated at the 
bottom right of each plot. The line has a slope of unity. Lower panel: as upper panel, but investigating the influence of using our 
radial velocities (rv MAX) or the SEGUE ones (rv SEGUE). Offset and scatters are calculated from the difference "SEGUE-MAX". 



4.1. Whole domain vs. spectral windows 

The classical x 1 fitting procedure uses every datapoint; there- 
fore, a straight-forward method to speed up the analysis would 
be to mask those parts of the spectrum which do not contain in- 
formation about the parameters, essentially those that are merely 
continuum. Certainly, by considering only the Call , Balmer and 
Mglb lines, the number of operations becomes smaller, increas- 
ing the processing speed ofx 2 calculations. This is clumsy how- 
ever: an empirical decision must be made about the relevance of 



pixels, and no extra weighting is considered - and the time taken 
for the parameters estimation is still long if one decides to do it 
for many spectra. The use of the b-vectors in the MA^ method 
means that very little weight is placed on pixels which do not 
significantly change with the parameter under study, automati- 
cally removing the sectors without lines. 

The remarkable result of the MA^ method is that it is possi- 
ble to determine the maximum of the compressed likelihood in 
10 milliseconds with a present day standard desktop PC. This is 
at least 300 times faster than the same procedure when doing an 
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efficient evaluation of the uncompressed data. 
The 1 cr confidence contours indicate errors in the parameters of 
0.24 dex in metallicity on average, 130 K in temperature and 0.5 
dex in gravity. Examples of fits are shown in Fig. [5] The upper 
plot shows the fit between a randomly selected SEGUE spec- 
trum and the best model - the crosses correspond to the observed 
spectrum and the dashed red line to the model, with stellar atmo- 
sphere parameters in the legend, as well as the resulting reduced 
X 2 of the fit. The plot in the bottom is the fit of the same star, 
but considering only the data points within the line regions iden- 
tified in Table Q] as discussed above. Again, the red dashed line 
indicates the model with parameters in the legend. The x 2 of the 
fit is also given. 

Because the lines contain the most information about metal- 
licity, temperature and gravity, we expect to obtain the same 
results whether we use all the data points or only those corre- 
sponding to the lines. The legends in Fig.[5]show the parameters 
estimated in both analyses. The small differences between them 
are within the 1 cr errors. We plotted in the upper panel of Fig. [6] 
the comparison for metallicity (left panel), temperature (middle) 
and gravity (right) of our sample of stars when using the whole 
spectral range ("whole") or only spectral windows with the lines 
("win"). The offsets and scatters of the distribution as well as 
the one-to-one relation are indicated in the legend of each plot. 
Here we randomly selected 150 stars to plot in Fig. [6] to visu- 
alize better the correlations of the results. Offsets and scatters 
are calculated with the entire star sample. Metallicity has excel- 
lent agreements, with a small scatter of 0.094 dex, as shown in 
the legend of the plot. Temperature has a negligible offset of 38 
K and a scatter of 69 K, which is also less than the accuracies 
obtained in the temperature estimation. Gravity has the largest 
scatter offset of 0. 1 dex and a scatter of 0.34 dex, but this is still 
within the errors. 

The negligible offsets obtained when using the entire spec- 
tral range against the limited regions is an encouraging result in 
terms of the effect of the G-Band on the spectrum. As pointed 
out above, the dependence on this feature in the parameters un- 
der study is not as strong as the rest of the lines of TableQ] This is 
translated into less weight for the compression, as seen in Fig.Q] 
meaning the G-band does not play a crucial role in our com- 
pressed data set for the parameter recovery. 

Let us remark that this does not mean a lack of CEMP stars in 
our sample, we simply do not see them in the compressed space. 
A spectrum with a strong observed CH molecule gives a similar 
compressed^- 2 as one with a weak one, because the compression 

does not consider variations in carbon abundance. The full x 2 , on 
the other hand, will certainly be larger for the spectrum with a 
strong G-band, because our models do not have different carbon 
abundances. 

In any case, the fraction of CEMP stars is high and it is cer- 
tainly interesting to locate them with the MAy method in the 
future. This implies that we need to create another dimension in 
the grid of synthetic spectra - models with varying carbon abun- 
dances - and increase the number of parameters to analyze to 
five. One more dimension certainly means a heavier grid of mod- 
els, but in terms of parameter recovery, there is no big difference 
to use four or five y-numbers for the compressed^ 2 calculation. 
This study goes beyond the scope of this paper, however, where 
we only introduce the MAy method. 



4.2. Effect on radial velocities 

As mentioned in Sect. 13.31 and seen in Fig. |4] radial velocities 
given by the SEGUE database have a mean difference of 20 km/s 
compared to ours. To study the effect of this difference in the fi- 
nal parameter estimation, we analyzed a sub-sample of randomly 
selected 2,670 stars of our sample using SEGUE radial veloc- 
ities. The comparison of the metallicity, effective temperature 
and surface gravity obtained when using our radial velocities 
("rv MAX") and the SEGUE ones ("rv SEGUE") can be seen 
in the lower panel of Fig. [6] In this plot, also a small sample 
of randomly selected 150 stars was used to better visualize the 
results, but offsets and scatters were calculated considering the 
2,670 stars. A difference in 20 km/s produces offsets and scatters 
in the three parameters, with 0.007 + 0.040 dex for metallicity, 
-5.08 ± 30.79 K for temperature and 0.01 + 16 dex for gravity. 
These offsets are negligible when compared with the lcr errors 
obtained the parameter estimation. 

4.3. Comparison with the SEGUE Stellar Parameter Pipeline 

The S E GUE Stellar Parameter Pipeline (SSPP, iLee et all 
2008a b; Allende Prieto et al. 2008) is a combination of different 
techniques to estimate the stellar parameters of SEGUE. Some 
of them are to fit the data to a gri d of synthetic spectra like 
the k24 (lAllende Prieto et. alJ I2006I) . and the kl3, NGS1 and 
NGS2 ones dLee et alj|2008al). Another m ethod are the ANNSR 
and ANNPvR of lRe Fiorentin et all d2007l) . which are an artificial 
neural network (ANN) trained to use a grid of synthetic (S) or 
real (R) spectra to estimate t he parameters of real spectra (R). 
Other options are line indices (fWilhelm et alJI 1999h and the Call 
K line index dBeers et alJll999h . 

We compare our results with the final adopted SSPP value 
(ADOP), and those of the k24, NGS 1 and ANNSR grids of syn- 
thetic spectra. Each of these results can be found in the SSPP 
tables of the SDSS database. The grids of models are created us- 
ing Kurucz stellar atmosphere models, but k24 and NGS 1 cover 
the wavelength range of [4400,5500] A. The ANNSR grid in- 
cludes the Call triplet close to 8600 A. 

Figure [7] shows the correlations between the results of these 
methods and our own. As in Fig. [6] we plotted a random selec- 
tion of 150 stars to better visualize the correlations, but values of 
offsets and scatters were made using the entire sample of 17,274 
stars. The left panels correspond to the metallicity, the middle 
ones to temperature and the right ones to the surface gravity. The 
top row of plots shows a general comparison of all the different 
methods, plotted with different colors. The y-axis indicates the 
ADOP results and the x-axis the methods MA^ (black), NGS1 
(blue), ANNSR (red) and k24 (green). The one-to-one line is 
also plotted in the figures. For each distribution a histogram of 
(methodj - methodj) with i + j was fitted with a Gaussian to 
obtain estimations for offset (mean A) and scatter (standard de- 
viation cr). The results are summarized in Table [2] In the lower 
panel of Fig. [7] we plotted the comparisons of individual meth- 
ods for 150 randomly selected stars. These panels help to better 
visualize how the single methods of the pipeline and our own 
correlate with each other. Figure|7]and Table|2]show that our re- 
sults reasonably agree with the pipeline, and that the scatter is 
of the same order on magnitude as that between the individual 
SSPP methods. 

Considering individual parameters, our metallicities have a 
general tendency to be ~ 0.3 dex lower than adopted metal- 
licity of SSPP (ADOP). A similarly large offset exists between 
ANNSR and NGS1, with -0.22 dex. This could be due to the 
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Fig. 7. Upper pl ot: Comparison o f the results fro m the SEGUE Stellar Pa rameter Pipeline (SSPP) with different method s: This work 
(MAX, black), dLee et al.l2008al NGS1, blue), dRe Fiorentin et al.ll2007l ANNSR, red) and dAllende Prieto et. alJ2006l k24, green). 
Lower part: Comparisons of individual methods for a randomly selected subsample of stars. Each plot has a line with slope of unity. 
The different offsets and scatters between the results are indicated in Table|2] In each plot, a selection of 150 random stars has been 
used. 



consideration of the Call lines in the fit, which are not included 
in the NGS1 and k24 ones. ANNSR also includes Call lines, the 
offset of -0.17 dex to our results is smaller than 1 <x error. If the 
same lines are used in the fitting, no offset should be seen, as 
will be shown in Sect. 



The temperature, on the other hand, shows a small offset of 
-61 K with respect to the ADOP value of the pipeline. It is en- 
couraging to see that we obtain the best agreement, except for 
ANNSR, which has no offset. A large mean difference is found 
between k24 and NGS 1 , where the offset is 244 K. The scatter in 
the various comparisons varies from 100 K (NGS1 v/s ADOP) 
to 200K (k24 v/s ANNSR); our scatter of 1 12 K with respect to 
ADOP is one of the lowest values. 



Finally, logg shows the largest offset. We derived gravity 
values 0.51 dex higher than the ADOP ones. The worst case is 
the comparison between our method and ANNSR with 0.63 dex. 
Between the methods of the pipeline, the largest difference (0.27 
dex) is found between ANNSR and NGS 1 . The best agreement 
we found for gravities are with NGS1, with an offset of 0.38 dex. 
The scatter for gravities varies from 0.23 dex (NGS1 v/s ADOP) 
to 0.48 dex (ANNSR v/s k24 and MA^v/s k24). A possible 
approach to correct our gravities would be to shift the zero point 
by -0.38 dex to agree with those of NGS1. We prefer to accept 
that gravity is our least constrained parameter, given the lack 
of sensitive features except for the wings of the Mglb lines, 
and the noise in the spectra does not restrict logg effectively. 
In order to deal with this problem, k24 and NGS 1 smooth the 
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spectra to half of the resolution, gaining signal-to-noise by this 
procedure. We remark that the k24 grid includes g-r colors, 
also to constrain the temperature. This automatically leads to 
different gravity values. A more extensive discussion of this 
aspect will be given in Sect. 15.61 The right panels of Fig. [7] 
shows that all log g estimates are rather uncorrected with each 
other with respect to the other ones. 



4.4. Systematic errors: Choice of fiducial model 

The b-vectors are basically derivatives with respect to the pa- 
rameters. With the assumption that the responses of the model to 
these derivatives are continuous, the choice of the fiducial model 
(FM) used to build the b-vectors is free. This assumption, how- 
ever, is only correct for certain regions of the parameter space. 
For example, Balmer lines are good te mperature indicators for 
stars no hotter than 8000K dGravll 1992b . and for cold stars spec- 
tral lines are affected by molecular lines. Limiting our candidate 
stars to the range of 5000K to 8000K we know that we have a 
clean spectra with Balmer lines as effective temperature indica- 
tors. This allows us to choose a FM with any temperature in that 
range to construct our weighting vector bj e ff, which will repre- 
sent well the dependence of the model spectrum on temperature. 
Similar assumptions are made for the other two model parame- 
ters. 

In order to test this assumption we studied the effect of using 
different fiducial models in the final results of the parameter 
estimation. In Fig. [8] we plot the results of metallicity (left), 
temperature (middle) and gravity (right) for 150 randomly 
selected SEGUE stars. There are four sets of different fiducial 
models [A,B] used for the b vector calculation. The parameters 
of the FM are indicated at the right side in each set. From top to 
bottom the first one compares the results of a cold (T e fr = 5500 
K) and a hot (T e ff = 7500 K) fiducial model of our grid of 
synthetic spectra, the second one a metal-poor ([Fe/H] = -2.0) 
and a metal-rich ([Fe/H] = -0.8) one. The third one considers 
two different gravities (logg = 3.7 and logg = 4.7). For 
these three sets, all other parameters are identically calculated. 
Finally, the last one uses models which vary all three parame- 
ters. The histograms of each plot correspond to the difference 
of the results obtained with the fiducial model with a set of 
the parameters A and B, respectively. The Gaussian fit of the 
histogram is overplotted and its mean (A) and standard deviation 
(cr) are indicated in the legend at the top. 

Metallicity shows a well defined behavior under different b- 
vectors. There are negligible offsets when varying only one pa- 
rameter in the FM, except for the case where the FM is totally 
different. But even in that case the offset of -0.15 dex is less 
than the 1 cr errors of 0.25 dex obtained in the metallicity es- 
timation. We conclude that this parameter does not show real 
systematic offsets due to the choice of fiducial model. The scat- 
ter obtained in the metallicity is also less than the 1 cr errors of 
the results. Only in the last case (bottom left) it becomes com- 
parable with the errors. Because gravity is a poorly constrained 
parameter (see discussions above) and the y-numbers are un- 
correlated, the effect when using different gravities in the FM 
should be completely negligible in the result of the other param- 
eters. This can be seen in the third panel from top to bottom of 
Fig. [8] where metallicity and temperature show extremely small 
offsets and scatters. 

The derivation of temperature using different fiducial mod- 
els shows a similar behavior. The discrepancies between results 



Results for 150 stars with fiducial models A— B 
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Fig. 8. Histograms of the differences of the results obtained using 
8 different fiducial models for the compression. Left: metallicity, 
middle: temperature and right: gravity. The parameters of the 
fiducial model are indicated in the right side of each plot, labeled 
as [[Fe/H], T e fr, logg]. Each histogram has a Gaussian fit, where 
the peak (A) and the standard deviation (cr) are indicated in the 
top legend. 



when using different FM are on the order of 50 K, except when 
the metallicity is varied. The mean difference of the final tem- 
perature estimates when using a metal-poor fiducial model and 
a metal-rich one is of 141 K, which is smaller than some dif- 
ferences seen in Table [2] between the methods of the SEGUE 
Stellar Parameter Pipeline. This probably happens because the 
spectrum of a metal-rich star presents more lines than a metal- 
poor one, which is translated into the b-vector as temperature- 
dependent regions that do not exist. In the case of SEGUE data, 
many of the weak lines seen in a metal-rich synthetic spectrum 
and(or) the b-vector of a metal-rich fiducial model are hidden by 
the noise in the observed spectra and the assumption of a tem- 
perature dependence due to these weak lines may not be correct. 
For the analysis of low signal-to-noise spectra, it is preferable to 
choose a FM with a rather low metallicity. But even when using 
a high metallicity for the FM, the offset is comparable with the 
1 cr error of 150 K for the temperature estimations. We conclude 
that temperature does not show a significant systematic offset 
due to choice of fiducial model either. 
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Table 2. Offset s (A) and scatter (<x) of the differences (raw - colum n) betwee n the methods: MAy (this work). ADOP dLee et al.l 
l2008ah . NGSl (iLee et al.ll2008al) . ANNSR dRe Fiorentin et al.ll2007l) and k24 (lAllende Prieto et. al. 2006) for the entire sample of 
selected SEGUE stars. 



[Fe/H] 


ADOP 

A/a 


NGSl 

A/a 


ANNSR 

A/a 


k24 

A/a 


MA^ 
k24 
ANNSR 
NGSl 


-0.32/0.23 
0.04/0.14 
-0.13/0.22 
0.08/0.17 


-0.41/0.29 
-0.07 / 0.25 
-0.22 / 0.23 


-0.17/0.33 
0.17/0.27 


-0.34 / 0.26 


T e ff 




MA^ 
k24 
ANNSR 
NGSl 


-61/112 
98/101 
1/114 
-142/93 


80/132 
244/ 166 
147/173 


-76/ 170 
88 / 195 


-159/160 


l°g£ 




MA^ 
k24 
ANNSR 
NGSl 


0.51/0.39 
0.08/0.35 
-0.14/0.36 
0.09 / 0.23 


0.38 / 0.43 
-0.02 / 0.42 
-0.27 / 0.45 


0.63 / 0.47 
0.26/0.48 


0.41/0.48 



Finally, gravity shows larger discrepancies in the final results 
when varying the FM. The differences can be as large as 0.4 
dex in the worst case when using different metallicities in the 
FM. This could also be because a metal-rich FM contains more 
lines and therefore greater sensitivities to gravity, which in our 
low signal-to-noise spectra is not the case. These sensitivities are 
confused with the noise in the observed spectrum. Scatters of 0.5 
dex are in general on the order of the lcr errors. 

5. Application to high-resolution spectra 

As a check, MAx a l so was tested with a smaller sample of 28 
high-resolution spectra. For consistency with the low resolution 
implementation, again metal-poor dwarf stars were selected. 

5.1. Data 

The spectra w ere obtained with U VES, the Ultraviolet Echelle 
Spectrograph dDekker et aU EOQO) at the ESO VLT 8 m Kueyen 
telescope in Chile. The resolving power is R ~ 43,000 and the 
signal-to-noise is typically above 300 in our spectra. Because 
MAy is an automatic fitting tool for synthetic spectra, we first 
had to be sure that the observed spectra could be fitted by our 
grid of synthetic spectra. This means that we had to have easily 
distinguishable unblended spectral lines and a clear continuum. 
For these reasons we selected a part of the red setting of the 
UVES spectra, which covers the wavelength range 580 nm - 680 
nm, where there are many unblended lines. 

5.2. Grid of high-resolution models 

When modeling spectra in high resolution, there are some dif- 
ferences to the low-resolution case and the following considera- 
tions have to be made: 

- A value for micro turbulence must be set, because the shape 
of strong lines sensitively depends on the value chosen for v t . 
This is not the case in low resolution, so that there is no need 
to fine-tune this parameter. In high resolution the set of basic 
stellar atmosphere parameters is ([Fe/H], T e ff, logg, v f ). To 



do a proper parameter estimation we would need to create 
a four-dimensional grid of synthetic spectra, varying all the 
parameters indicated above. This goes further than the pur- 
poses of this paper, where we aim to check the applicability 
of our method. Therefore we fixed the microturbulence pa- 
rameter to a typical value, guided by the results of a standard 
"classical" spectral analysis (see below). 

- To compute synthetic spectra, a list of lines with their wave- 
length and atomic data must be included. The atomic data 
consist of oscillator strength log gf and excitation poten- 
tials, which are an important source of differences in the final 
shape of a given line. For the comparison consistent line lists 
are needed. 

- In order to make a proper fit to a broad wavelength range, the 
models should be perfect. Our models were computed under 
the simplifying assumption of LTE and with plane-parallel 
atmosphere layers; the atomic data have errors, too. There 
are too many features to fit, therefore a global fit becomes 
difficult, almost impossible. For this reason we selected win- 
dows of a limited spectral range and fitted the observed spec- 
tra only in these. 

We determined the three parameters temperature, metallicity 
and gravity using neutral and ionized iron lines. Neutral lines 
are sensitive to tem perature and single ionized ones to grav- 
ity (jFuhrmann 1998). The metallicity was effectively obtained 
from both Fel and II, for a given temperature and gravity. In the 
wavelength range of 580nm - 680nm, we have six Fell lines 
and 20 Fel ones that are unblended and strong enough to be 
present in most metal-poor stars of our sample. The lines are 
indicated in Table [3] where wavelength and log gf values are 
listed. The log g f v alues and excitation pote ntials are taken from 
iNissen et all (120021) and the VALD database dKupkaet al .12 000). 

Our grid of high resolution models covers a range in metal- 
licity of -2.5 < [Fe/H] < -0.5, in effective temperature of 
5500 < T e ff < 6600, in surface gravity of 3.5 < logg < 5 and in 
scaling factor for normalization from 0.85 to 1.15. The grid steps 
are the same as in Sect. 13.21 The models have [a/Fe] = +0.4 and 
v, — 1 .2 km/s, according to the values obtained for our stars with 
the "classical" method (see below). 
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Fig. 9. Fit of the spectrum of the star HD195633 (points) with the best model (red dashed line). The parameters of the synthetic 
spectrum are [Fe/H] = -0.5, T e a = 6318 K and logg = 4.61, with^- 2 = 6.31. The blue line corresponds to the model with [Fe/H] = 
-0.6, T e ff = 6005 K and logg = 3.86 (x 1 = 7.31), obtained with the "classical" approach. Each panel represents a spectral window 
where the lines from Table [3] are labeled at the bottom. 



5.3. Preparing the observed data for synthetic spectral fitting 



The reduction of the data was carried out in the same way 
as in Sect. 13.31 with the following difference in the normal- 
ization procedure. At this resolution, neither the polynomial- 



fitting approach bv lAllende Prieto et. al .1 (120061, norlRArf 5 ! were 
able to automatically find a good pseudo-continuum. Hence, 



IRAF is distributed by the National Optical Observatory, which is 
operated by the Association of Universities of Research in Astronomy, 
Inc., under contract with the National Science Foundation. 
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Table 3. Fel and II lines used for the high-resolution fits. The 
first column indicates the wavelength; the log gf value is taken 
from the VALD database and is given in the second column. 



wavelength (A) 


log gf 


Fel 


6136.62 


-1.410 


6137.69 


-1.403 


6191.56 


-1.416 


6219.28 


-2.448 


6230.72 


-1.281 


6246.32 


-0.877 


6252.56 


-1.767 


6265.13 


-2.540 


6335.33 


-2.177 


6336.82 


-0.856 


6393.60 


-1.576 


6411.65 


-0.717 


6430.85 


-1.946 


6494.98 


-1.239 


6592.91 


-1.473 


6677.99 


-1.418 


Fell 


6084.11 


-3.881 


6149.26 


-2.724 


6238.39 


-2.754 


6247.56 


-2.329 


6416.92 


-2.877 


6432.68 


-3.501 



the observed spectra w ere normalized interactively using Midas 
dCrane & B anse 1982), and the synthetic spectra were computed 
with normalized flux. The compression procedure was the same 
as in Sect. 13.41 Compressing the grid of high-resolution spectra 
takes more time than for the low-resolution ones, but this has to 
be carried out only once to obtain a new grid of y-numbers. 

5.4. Results of the high-resolution analysis using MAx 

Figure [9] shows an example fit of HD195633 (points) with the 
best MAjf model (red dashed line) and the model with param- 
eters found with the "classical" method (blue dashed line; see 
below). Each panel represents a spectral window with the lines 
of Table [3] used for the parameter estimation. All three parame- 
ters were determined simultaneously and the best fit corresponds 
to the model with parameters of [Fe/H] = -0.5, T e fr = 6318 K and 
logg = 4.61. The final parameters of all stars in this sample are 
given in Table|4] The first column of the table indicates the name 
of the star and the six next ones are the parameters found with 
two different MA^ analyses. The first set ("free") corresponds to 
the standard approach of determining all parameters simultane- 
ously from the spectrum only. In the second variant ("restricted") 
the gravity is determined independently using Eq. (14), which 
will be introduced below. The last four columns are the parame- 
ters obtained from the "classical" approach (below). 

5.5. Results of a "classical" analysis 

In order to compare our MAx results, we analyzed the spectra 
with a classical procedure, determining the stellar parameters 
through a n iterative process. The method is basically the same 
as in, e.g.. INissen et al.l (|2002h : 

For the effective temperature it relies on the infrared flux 
method, which provides the coefficients to convert colors to ef- 



fective tejnrjeratirres. Weused the (V-K) color and the calibra- 
tion oflAlonso et all (1 19961) . after converting t he K 2MASS filter 
(Skrutski e'et all 1 19971) to the Johnson filter dBessellll2005l) and 
de-reddeni ng the color. The ex tinction was taken from Schlegel's 
dust maps (Schl egel et al.ll 998) in the few cases where we could 
not find the values in Niss en et al. I d2002ll200l . 

Surface gravities were determined using the basic parallax 
relation 

logg = log — + log + 0.4 (Vc + BC) 

Mq 1 



+2k>g7T + 0.12 + log g e , 



(14) 



where Vc is the V magnitude corrected for interstellar absorp- 
tion, BC the bolometric correction and n the parallax in arcsec. 
W e adopted a differen t mass value for each star based on those 
of iNissen et al.l d2002l) . which are between 0.7 and 1.1 M . The 
bolometric correction B C was calcula t ed wi th the solar calibra- 
tion of BC Q = -0. 12 as in INissen et all d 19971) . 

The equivalent widths of neutral and ionized iron lines were 
used to de termine the metallici ty [Fe/H]. It was obtained using 
Fitline dFrancois et alj2007h . which fits Gaussians to the line 
profiles. The equivalent widths computed from these G aussians 
fits w ere converted to abundances by running MOOG dSnedenl 
fl97l Sobeck priv. comm. 200 8) with the plane-paral lel LTE 
MARCS model atmospheres of dGustafsson et al.ll2008l) . which 
differ from those used for MAjf . The line list contains the lines 
of Table [3] and more taken from bluer wavelengths in the range 
300 - 580 nm (Hansen, priv. comm. 2009). 

Finally, the value for the microturbulence was found by re- 
quiring that all the equivalent widths of neutral lines should give 
the same Fe abundance as the ionized lines. 

In order to obtain the four parameters in this way, we started 
with an initial guess for each of the interdependent parame- 
ters. After determining them with the above steps all values will 
change due to the interdependence, hence we had to iteratively 
determine the parameters until their values showed a negligible 
change. 

We are aware that our metallic ity based on Fel lines is i gnor- 
ing any non-LTE effects dAsplundll2005l; ICollet et alj|2005l) . But 
because neither the classically nor the MA^ analysis are taking 
this into account, the results obtained by both methods can well 
be compared. 

The upper panels of Fig. \W\ show the results of metallic- 
ity (left), temperature (middle) and gravity (right) for the entire 
sample of stars. The jc-axis shows the results of the automatic 
fitting MAx an d the y-axis the "classical" parameter estimation 
(EW). The one-to-one line is overplotted in each figure and the 
legend indicates the mean of the differences (MA^ - EW) and 
its standard deviation, denoted as offset and scatter, respectively. 
Each point corresponds to a star of Table [4] with parameters ob- 
tained by the "free" method. 

Gravity shows a scatter of <x = 0.46 dex and a negligible off- 
set. The determination o f the gravity from Fell lines has always 
been a problematic task: Fuhrmann] d 19981) in an extensive study 
of parameters of nearby stars showed that surface gravities of F- 
type stars located at the turn-off point can easily differ by up to 
~ 0.4 dex if derived from either LTE iron ionization equilibrium 
or parallaxes. The amount of ionized lines (gravity dependent) 
present in the spectra is usually smaller than that of neutral ones 
(temperature dependent). By performing an automatic fitting of 
weighted spectra that contain only six Fell lines compared to 20 
Fel lines, a scatter of 0.46 dex is reasonable. 

The result for the effective temperature shows a very small 
offset of only -10 K, but a quite large scatter of <x = 227 K. 
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Table 4. Parameters of the stars (Col 1) obtained with MA^ for both types of analysis discussed in the text. In the first analysis 
all three parameters are determined with only the spectrum, while in the second one log g was determined independently by using 
parallaxes according to Eq. (14), as for the classical method (EW; last four columns) 





MAx, free 




MA^, restr. 






classical, EW 




star 


Teff (K) 




[Fe/H] 


7W (K) 




[Fe/H] 


Teff (K) 




[Fe/H] 


v, (km/s) 


HD3567 


6103.5 


4.79 


-1.200 


6176.0 


4.20 


-1.223 


6035.0 


4.08 


-1.200 


1.50 


HD 19445 


_ 


_ 





5746.5 


4.10 


-2.293 


5982.0 


4.38 


-2.075 


1.40 


HD22879 


5770.6 


4.31 


-0.877 


5778.1 


4.31 


-0.880 


5792.0 


4.29 


-0.885 


1.20 


HD25704 


5638.9 


3.61 


-1.094 


5864.4 


4.19 


-0.912 


5700.0 


4.18 


-1.010 


1.00 


HD63077 


5795.5 


4.21 


-0.796 


5787.7 


4.21 


-0.793 


5629.0 


4.15 


-0.935 


0.90 


HD63598 


6078.6 


4.79 


-0.681 


5804.1 


4.23 


-0.901 


5680.0 


4.16 


-0.885 


0.90 


HD97320 


5601.1 


4.65 


-1.411 


6025.1 


4.24 


-1.183 


6105.0 


4.18 


-1.175 


1.20 


HD 103723 


5999.2 


4.29 


-0.896 


6004.4 


4.29 


-0.896 


6128.0 


4.28 


-0.805 


1.50 


HD 105004 


5802.5 


4.33 


-0.900 


5802.5 


4.33 


-0.900 


5900.0 


4.30 


-0.800 


1.10 


HD106038 


5994.6 


4.48 


-1.392 


5995.5 


4.48 


-1.392 


5950.0 


4.33 


-1.390 


1.10 


HD 113679 


6037.3 


4.73 


-0.500 


5757.5 


4.35 


-0.708 


5759.0 


4.04 


-0.590 


0.90 


HD 11 6064 


5757.7 


4.58 


-2.41 


5757.7 


4.58 


-2.401 


5999.0 


4.33 


-2.000 


1.50 


HD121004 


5509.2 


4.05 


-0.670 


5657.8 


4.40 


-0.650 


5711.0 


4.46 


-0.725 


0.65 


HD122196 


5811.8 


3.50 


-1.954 


5810.2 


3.56 


-1.954 


6048.0 


3.89 


-1.760 


1.20 


HD 126681 


5556.9 


4.44 


-1.055 


5521.6 


4.68 


-1.075 


5532.0 


4.58 


-1.170 


0.60 


HD160617 


5804.1 


3.50 


-1.961 


5985.2 


3.75 


-1.794 


6028.0 


3.79 


-1.785 


1.30 


HD175179 


6080.5 


5.00 


-0.500 


5779.5 


4.47 


-0.691 


5758.0 


4.16 


-0.690 


0.90 


HD188510 


5568.8 


4.35 


-1.341 


5500.0 


4.65 


-1.231 


5536.0 


4.63 


-1.600 


1.00 


HD 189558 


5610.8 


3.58 


-1.191 


5500.0 


3.75 


-1.181 


5712.0 


3.79 


-1.145 


1.20 


HD195633 


6318.2 


4.61 


-0.500 


6037.7 


3.99 


-0.687 


6005.0 


3.86 


-0.625 


1.40 


HD205650 


5693.7 


4.93 


-1.195 


5754.0 


4.47 


-1.208 


5842.0 


4.49 


-1.170 


0.90 


HD298986 


5891.5 


3.66 


-1.592 


6297.9 


4.22 


-1.304 


6144.0 


4.18 


-1.410 


1.40 


CD-333337 


5808.2 


3.66 


-1.598 


5946.9 


3.85 


-1.490 


5952.0 


3.95 


-1.435 


1.40 


CD-453283 


5500.0 


4.26 


-0.978 


5885.5 


5.00 


-0.685 


5657.0 


4.97 


-0.835 


0.80 


CD-571633 


5939.9 


4.37 


-0.893 


5939.9 


4.37 


-0.893 


5907.0 


4.26 


-0.930 


1.10 


CD-3018140 


6143.1 


3.50 


-2.095 


6246.5 


3.50 


-1.980 


6340.0 


4.13 


-1.870 


1.00 


G005-040 


6228.2 


5.00 


-0.584 


5930.4 


4.52 


-0.794 


5766.0 


4.23 


-0.840 


0.80 


GO 13-009 


6529.5 


4.58 


-2.29 


6556.0 


4.08 


-2.299 


6416.0 


3.95 


-2.225 


1.40 



Given that two different methods were employed to determine 
it, this still appears to be acceptable. The scatter may also be af- 
fected by our fixed value for the microturbulence, which was 
taken from the average of the values found with the EW ap- 
proach (see Table but which in individual cases may differ 
severely. 

One can decrease the scatter in the temperature difference 
by removing the three most discrepant objects from the sam- 
ple. The star HD63598 (triangle) has an offset of 398 K. The 
Schlegel maps give a de-reddening for this star that is unrealisti- 
cally large, so ins tead we set it to zero . This led to a temperature 
that was too low. iJonsell et al.l d2005l) have found a temperature 
of 5845 K for this star, reducing the difference to 233 K with 
respect to our result. The star G005-040 (asterisk) is the sec- 
ond case, where the offset is 462 K. For this star, the continuum 
subtraction was not perfect in every spectral window, resulting 
in a fit where the parameters were quite unreliable. Finally, the 
weak lines of the observed spectrum of HD 19445 pushed MA^ 
to the border of the grid. The parameters in this case were un- 
determined. By removing these three stars from our sample, the 
scatter for temperature is reduced to 197 K. 

Metallicity, on the other hand, shows a very good agreement, 
with a negligible offset of 0.02 dex and a scatter of 0.16 dex 
(upper left panel of Fig. HOI . The lines used for the automatic 
fitting (MA^) and for the classical analysis (EW) are in most 
cases identical, except for the lowest metallicities, where some 
of them are hardly visible. In these cases the classical method 



also resorted to other lines outside the MA^ wavelength range. 
The atomic data are identical, but the value for microturbulence 
are not, as mentioned before. In view of all this, the very good 
correlation of metallicity is encouraging. 

5.6. "Restricted" parameter recovery 

An example for the individual fit of the best MA^ model to the 
observed spectra was shown in Fig. [9] for HD195633, which 
is shown as a filled symbol in Fig. [TpJ and is one of the ob- 
jects with the most pronounced discrepancies (see also Table©. 
Nevertheless, the two synthetic spectra appear to fit equally well 
the observations (x 2 = 6.31 for the MA^ model and^ 2 = 7.31 
for that obtained from the classical analysis). We could not de- 
cide from the x 2 values which method leads to a more accurate 
parameter determination for this star. One would need indepen- 
dent information of the spectrum about this star to reach conclu- 
sions about its parameters. 

It is also interesting to notice from the upper panels of Fig.fTOl 
that stars with large discrepancies in gravity (for example those 
with special symbols) also give large discrepancies in tempera- 
ture, as expected. They show a quite unsatisfying fit for the Fell 
lines. Motivated by this, we did an additional test by restricting 
MA^ to the determination of the parameters with input values 
for logg obtained from the EW method, i.e. we found a local 
maximum point of the likelihood in a restricted area. To do this 
we chose the three closest gravity values from our grid of syn- 
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thetic spectra to the classical one found with Eq. (TBI - see Col 9 
of Table|4]- and we searched for the maximum point within this 
range. May be that for log g there is no local maximum in this 
range, and the final estimation will go to the border of the grid, 
as the case of CD-3018140, where \ogg = 3.5. The general ten- 
dency is anyway a local maximum close to the input EW-value. 

Now the agreement with the classical method for the effec- 
tive temperatures became excellent, with a negligible offset of 
only 7.42 K and a scatter of 128 K, as seen in the lower pan- 
els of Fig. [10] Gravity was also better constrained, with a small 
scatter of 0.14 dex and a negligible offset. Note that we did not 
necessarily obtain identical values, because the MAx gravity is 
obtained by using the value from Eq. (TBI) only as input, and we 
were looking for a final solution close to this value. The behav- 
ior of the metallicity does not change with respect to the "free" 
case, demonstrating the robustness of the determination of this 
parameter. The results obtained when using parallaxes for the 
initial guess for logg are summarized in Table [4] in the columns 
under the heading "MAy, restr.". 

It is instructive to discuss the implication of this comparison. 
For the "classical" EW method we determined the parameters 
making use of the best information available and of the free- 
dom to adopt the method to each star individually. The iterative 
process allowed us to decide where to stop the iteration, or, if 
no satisfying convergence could be reached, to draw on the op- 
tions to move to another spectral window, to use other lines, or 
disregard problematic lines. Moreover, the continuum could be 
separately subtracted for each line, thus creating locally perfect 
normalized flux levels. For the effective temperature, the more 
reliable infrared flux method employing photometric data could 
be used instead of only relying on spectra. The advantage of de- 
termining log g independently was already demonstrated. 

On the other hand we were attempting an estimate of the pa- 
rameters only from the spectral information without fine-tuning 
the models or fit procedure for each star for our automatic MAj 
method,. Given this, the comparison with the full interactive 
method is surprisingly good. We demonstrated that by using ad- 
ditional information for one parameter (here logg) it becomes 
even better for the remaining two quantities. Table [4] shows our 
final parameter estimate for our stars, when we make use of the 
parallaxes as additional information. 

The result of this test is that we demonstrated our ability to 
estimate the basic stellar atmosphere parameters quickly and ac- 
curately enough for a substantial sample of stars. While individ- 
ually severe outliers may occur, the method is accurate overall 
for a whole population of stars. For this purpose the EW method 
would be much too slow and tedious. Our method may also serve 
as a quick and rough estimate for a more detailed follow-up anal- 
ysis. 



6. Summary and conclusions 

We described MAy, a new derivative of MOPEDfor the estima- 
tion of parameters from stellar spectra. In this case the parame- 
ters were the metallicity, effective temperature and surface grav- 
ity. The method reduces the data to a compressed data set of three 
numbers, one per parameter. Assuming that the noise is inde- 
pendent of the parameters, the compressed data contain as much 
information about the parameters of the spectrum as the entire 
data itself. As a result, the likelihood surface around the peak 
is locally identical for the entire and compressed data sets, and 
the compression is 'lossless'. This massive compression, with 
the degree of compression given by the ratio of the size of the 



data to the number of parameters, allows the cost function cal- 
culation for a parameter set to be sped up by the same factor. 
For SDSS data, the spectra are of the order of 1000 datapoints, 
each with a corresponding error, which must be compared to a 
similarly sized model. The speed up factor in this case is then at 
least 1000/3, ~ 333x. 

This extremely fast multiple parameter estimate make the 
MA^ method a powerful tool for the analysis of large samples 
of stellar spectra. We have applied it to a sample of 17,274 metal- 
poor dwarf stars with low-resolution spectra from SEGUE, us- 
ing a grid of synthetic spectra with the parameter range of 
[-2.5, -0.5] dex in metallicity, [5000,8000] K in effective tem- 
perature and [3.5, 5] dex in surface gravity, covering a wave- 
length range of [3850, 5200] A. 

From the Call , Balmer and Mglb lines, which are the 
strongest absorption features identified in SDSS spectra, we es- 
timated the metallicity with averaged accuracies of 0.24 dex, the 
temperature with 130 K and logg with 0.5 dex, corresponding 
to the lcr errors. Surface gravity is a poorly constrained parame- 
ter using these data, mainly due to the lack of sensitive features 
of this parameter (apart from some degree of sensitivity of the 
wings of the Mglb triplet) and the considerable noise in the 
spectra when compared to high-resolution spectra. Additional 
information to the spectra, such as photometry, would help to 
constrain the gravity parameter more. 

MA^ has the option to simultaneously fit different spectral 
windows. We have compared estimates of the parameters using 
the whole spectrum and only those data ranges where known 
lines exist. Both analyses take approximately the same time and 
agree excellent in recovered parameters. This suggests that for 
these low-resolution spectra there is no need to carefully and la- 
boriously select the spectral windows to be analyzed with MAf, 
the method calculates this automatically as part of its weighting 
procedure. 

With the assumption that the spectra behave similarly under 
changes of the parameters, the choice of the fiducial model for 
the compression is free. We have created eight compressed grids 
using different fiducial models for the b-vector calculation and 
obtain agreement between the recovered parameters. Given the 
low signal-to-noise of our data, it is better to use a more metal- 
poor fiducial model, as the dependence on the parameters will 
be focussed on the strong lines. 

We have comprehensively investigated the correlations of 
our result s with those obt ained for the SEGUE Stellar Parameter 
Pipeline dLee et al.ll2008ah . which reports results from a num- 
ber of different methods. The results from MA^ agree well with 
those of the various pipelines, and any differences are consis- 
tent with those between the various accepted approa ches. We 
are aware of the MATISSE dRecio-Blanco et al .1120061) method 
of parameter estimate, which uses a different combination of 
weighting data but is closer to the MA^ approach than the stan- 
dard pipeline methods. We look forward to comparing our re- 
sults with those of MATISSE when they become available. 

More specifically, temperature agree excellently with the av- 
eraged temperature of SSPP, with a negligible offset of -61 K and 
low scatter of 1 12 K. Our metallicities show a tendency to be - 
0.32 dex lower than SSPP averages. The small scatter of 0.23 dex 
suggests that the different spectral features used in the analysis 
(mainly Call lines) could shift the zero point of the metallicity. 
The most pronounced discrepancy is found in surface gravity, 
where Ma^ reports values 0.51 dex (cr = 0.39) higher than the 
averaged value of the pipeline. We saw that this is consistent 
with the discrepancies seen between other methods. 
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[Fe/H] MAX Teff MAX (K) logg MAX 



Fig. 10. Upper panels: Correlations between the results obtained with the automatic fitting (MAX) and the classical method (EW) for 
metallicity (left), temperature (middle) and gravity (right). In each panel, the one-to-one line is overplotted and a legend containing 
the mean difference (in the sense "EW-MAy") and its standard deviation. Here the three parameters were determined simultaneously 
and directly from the spectra by our method. The stars with special symbols are indicated in the legend and correspond to special 
cases where the offset between the two methods is particularly large (discussed in text). Lower panels: As upper panels, but MA^ 
results were obtained forcing the method to find the best-fitting model close to the gravities computed from the parallax formula, 
Eq. (14). Individual data can be found in Table [4] (columns "MAy, free" for the upper case and "MA^, restr." for the lower one). 



We have also tested MA^ on a sample of 28 high resolution 
spectra from VLT-UVES, where no offset in metallicity is seen. 
In this case we have carefully chosen the models and spectral 
range for comparing our parameter estimation against a "clas- 
sical" approach: temperature from photometry, surface gravity 
from parallaxes and metallicity from equivalent widths of neutral 
and ionized iron lines. We have calculated the parameters our- 
selves to avoid additional systematic offsets that various differ- 
ent methods would introduce in stellar parameter scales. These 
results were compared with the automatic fitting of 20 Fel and 6 
Fell lines (that coincide in most cases with the equivalent widths 
calculations) made with MA^. We obtained large scatter in grav- 
ity and temperature (0.44 dex and 220 K, respectively), but no 
systematic offset. The normalization of the continuum in some 
cases was not perfectly done, making it difficult to fit every line 
properly: especially for the Fell lines. They produce a scatter in 
gravity, which drives a scatter in temperature. 
We have seen in our fits that our best model does not differ very 



much from the model with parameters found by the "classical" 
approach. Motivated by this we decreased the scatter in grav- 
ity by using additional information to the spectrum, i.e. using 
the gravities determined from parallaxes as input value. This 
forced our method to find a local maximum point of the like- 
lihood close to this input gravity value. Fixing the gravity gave 
an improved agreement in temperature, now with a scatter of 128 
K. Metallicity does not change when forcing gravities, illustrat- 
ing the robustness of determination for this parameter. 

MAx is an extremely rapid fitting technique and as such is 
independent of model and data used. It will work for any star for 
which an appropriate grid of synthetic spectra can be calculated. 
Although the grid calculation is time consuming, it only needs 
to be performed once allowing an extremely rapid processing of 
individual stars. 

MA^ is one of the fastest approaches to parameter deter- 
mination, and its accuracy is comparable with other methods. 
To develop MAy further in preparation for the next generation 
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of surveys (eg. Gaia, APOGEE or LAMOST), we plan to ex- 
pand our range of model grid and integrate photometric data 
to improve the determination of the parameters, especially sur- 
face gravity. We must be aware that in a grid of synthetic spec- 
tra with a broader parameter range a set of b-vectors from one 
single fiducial model will not represent well the dependence on 
the parameters in the entire sample. To cover the entire sample 
with well represented dependences, we plan to compute differ- 
ent compressed grids, each of them with b-vectors calculated 
from different fiducial models. The parameter estimate in this 
case will be made with the different compressed grids using a 
final convergence test. After convergence, the final parameters 
should be estimated from the compressed grid with the fiducial 
model with parameters close to the final result. 
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